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Abstract —This letter proposes an analytical approach to for¬ 
mulate the power system oscillation frequency under a large 
disturbance. A fact is revealed that the oscillation frequency is 
only the function of the oscillation amplitude when the system’s 
model and operating condition are fixed. Case studies also show 
that this function is damping-insensitive and could be applied to 
an inter-area model of a multi-machine power system. 

Index Terms —Oscillation frequency, power system, large dis¬ 
turbances, pendulum system, damping, small signal analysis. 

I. Introduction 

P OWER systems’ electromechanical oscillations are detri¬ 
mental to system stability. A traditional small signal 
analysis approach linearizes a system’s nonlinear differential- 
algebraic equation model near a specific operating condition 
to estimate its natural oscillation frequency (OF) but neglect 
the nonlinearity of the system, so the natural oscillation 
frequency cannot describe the accurate OF, especially under 
a large disturbance. When estimating the OF from power 
system measurements, most signal processing based methods, 
e.g. Fourier Analysis and Prony Analysis, need data over a 
time window much longer than the period of the targeted 
oscillation mode to estimate an average OF while ignoring 
any instantaneous change in the window. For a pendulum 
system, it is feasible to analytically formulate its period or 
frequency under a large disturbance by elliptic integral of the 
first kindJT] but there is no similar study for power systems, 
even for a single-machine-infinite-bus (SMIB) or a two-area 
power system, which can be regarded as a fictitious damped 
pendulum system with a constant torque. 

This letter will analytically formulate the power system’s 
OF under a large disturbance and introduce the Frequency- 
Amplitude curve as a characteristic in nature with the power 
system about a specific oscillation mode. The formulation will 
then be tested on a SMIB system and a two-area system. 

II. Oscillation Frequency Formulation 
Considering the motion equation 

f AS = uj 0 Auj r 

\ AA r = (P m - P e - DAuj r )/2H ' 

wo is the synchronous frequency, AS is the rotor angle 
deviation relative to its steady-state value <5o, A uj r is the 
per-unit rotor speed deviation, P m = P max sin So and P e = 
Pmax sin(<5o+A<5) represent the per-unit mechanic and electric 
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torques, respectively, where P max is the steady-state maximum 
power transfer, and H and 1) represent the inertia and damping 
factor of the machine, respectively. For simplicity, assume P m 
is constant under the disturbance. Eliminating A uj r leads to 
©, where /? = P max uJo/2,H for simplicity. 

A<5 + ( D/2H)AS + /3 sin(<5o + AS) = /3 sincio (2) 

Unlike a pendulum system whose oscillation is symmet¬ 
ric about the equilibrium, a SMIB system has no such 
symmetry0: for each cycle of its oscillation, the motions in 
the upper and lower halves of cycle are different and hence 
should be formulated separately. Assume that the maximum 
and minimum angle deviations, A S ma x and AS rrlln , are known 
and the damping is zero. 

The upper half of cycle is the time taken from the equi¬ 
librium to the maximum plus the returning time. Only the 
upswing is considered according to the symmetry of oscil¬ 
lation within this half cycle. Define potential energy at the 
equilibrium to be zero. The kinetic and potential energies are 
calculated by ( 0 . From the law of conservation of energy, the 
total energy holds from Ac) to A S max during the upswing. 
Solve for df and get 0. 

f E k = {dA5/dt) 2 H/uo (3) 

{ E p = P max (cos <5 0 - cos(5 0 + A<5) - AS sin <5 0 ) 


At = [2/l(cos(A(i + 5 0 ) ~ cos(A<S max ^o) 

+ (A6 - ASmax) sin do)] 2 dA<5 (4) 

By integrating both sides for dA<5 from 0 to AS max , the left 
side is the total time of that motion, doubling which gives the 
time of the upper half of cycle 0 . Next, change the integral 
variable by 0 , where k = sin(A<5 mox /2). Then we have 0 . 


nA5m ax 

T u = 2 [2/3(cos(A6 + (5 0 ) — cos(A S max 

JO 

+<5o) + (A(5 - ASmax) sin(5 0 )] 2 dAd (5) 
sinvj = sin(A(5/2)/sin(A6 maa; /2) (6) 


AS = 2arcsin(/csinyj) 
dA(5 = 2fccosy>(l — k 2 sin 2 ip)~^Aip 

Note that when AS increases from 0 to AS max , ip increases 
from 0 to 7t/2. Thus, 0 becomes 

f K / 2 2k cos ip 

T u = 2 g u (simp) = Atp (8) 


\/l - k 2 


sin 2 ip 
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g u = [2/3(cos(2arcsin(/csiny>) + <5 0 ) — cos(A S max 

+<S 0 ) + (2 arcsin(fc sin <^) - A<5 max ) sin<5 0 )] 5 (9) 

Similarly, T) will be 

[~' K / 2 2k cos 02 

Ti = 2/ -ffl(smy) . =d(jg (10) 

Jo v 1 — fc 2 sin 

where gi has the same form as g u except that A 5 max is 
replaced by A S m in and k = — sin(A5 m j n /2). Then, define 
OF / as ([ID. 

f = l/(T u + T t ) (11) 


When AS m ax and AS rnrn . approach 0, the general OF 
defined by ( ITTb will degrade to the nature frequency of the 
linearized model of the power system. Here is 

lim / = V /3 cos 8q/2t: (12) 

ASmax ,A5min —*0 

The following two methods can either solve ASmax and 
A S m in directly from the model or obtain their values from 
measurement data: 

Assume initial values of (0 to be A<j(0) and A<j(0). 
Because = 0 at each extreme point A<5 ep (either A5 max 
or A Smin), the conservation of energy gives (fl3l > about A S ep . 
When the system is stable. (Qj} has one positive root and one 
negative root, i.e. ASmax and A<5 min , respectively. 

(Aj(0)) 2 + /3( cos(A S ep + S 0 ) - cos(Ad>(0) + S 0 ) 

+(A<5 ep — AJ(0)) sinb 0 ) = 0 (13) 


Alternatively, if the time series of <5 are measured, by 
subtracting an estimate of its steady-state value, we may obtain 
values of A S max and A S m in- 

Finally, two approaches can be employed to calculate the OF 
in CD. The first approach uses power series to approximate 
® and ( I I Oi l which gives the OF by the sum of the first N 
terms of the series as shown in d, where f u and fi are 
integrands of ([Sj and (ITOt . respectively. When only keeping 
the first term, the approximation of OF is shown in dT5l >. The 
second approach is to estimate (0 and (flOt numerically. 
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fps — 1/ 


f-x/2 


fu ] (0) sin* tp + f[ l) (0) sin* (-p) 


dp 


fi 


ps\N=l 


1=0 17 u 

VP /sin(A<5„ 


C /2) sin (—AS 77 


V2tt V m(AS max ) 


m(A Smin) 


(14) 
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m(x) = V cos — cos(5q + x) — x sin So 
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06 ) 


Define the oscillation amplitude (OA) as 


( A5max — A5rr 


Equation (fl4b discovers that OF of a specific mode is only 
the function of OA when the systems model and operating 
condition are fixed. Thus, a Frequency-Amplitude (F-A) curve 
can be drawn from (IT4l) as a characteristic of the system. The 
OA when OF approaches zero indicates the stability limit. 


III. Case Study 

First, the influence of damping is illustrated on a SMIB 
power system, and then the proposed general OF is validated 
on Kundur’s two-area power system®. 


A. Test on a SMIB Power System 

Let H = 3, wo = 1207T, P ma x = 1-3 and Sq = 0.8 in 0. 
Consider three cases: 1) D = 1, At5(0) = 30° and Aj(0) = 
2rad/s; 2) D = 1, A6(0) = 60° and AS{0) = 2rad/s; 3) 
D = 3, A<5(0) = 60° and Ai5(0) = 2rad/s. In each case, 
solve the trajectory of Ac) using the Runge-Kutta method, and 
for any adjacent two extreme points on the trajectory, their 
differences in time and angle respectively give estimates for 
the OF and OA. The formulation (ITTI) is calculated numerically 
by Simpson’s rule. 

Figd gives the F-A curves about 3 cases and the theoretical 
F-A curve from the formulation, which all curves match well. 
Comparison between cases 1 and 2 indicates that the initial 
value does not change the F-A curve. Comparison between 
cases 2, 3 and the formulation shows that the F-A curve is 
insensitive to damping. 

B. Test on the Kundur’s System 

Only the inter-area mode is considered and a SMIB equiv¬ 
alent of the system is considered by defining the equivalent 
angle difference as O- Case 1 adds a three-phase fault on bus 
7 and clears it after 10 cycles by tripping three lines between 
buses 7 and 8. Case 2 extends the fault duration to 19.9 cycles 
to approach the transient stability limit. Fig[2]indicates that the 
proposed formulation could also be applied to a multi-machine 
power system. 

. r S\H\ + S 2 H 2 SoHo+SiHn 

A 0 = -(1/) 

Hi+H 2 fTs + fti 



Oscillation Amplitude/deg 



Fig. 1. Results on a SMIB 


Fig. 2. Results on Kundur’s system 


IV. Conclusion 

Analytically formulation for oscillation frequency of a 
power system has been proposed as a function of the oscilla¬ 
tion amplitude when the system’s model and operating condi¬ 
tion are fixed. The Frequency-Amplitude curve is introduced 
to characterize the system’s nonlinearity with an oscillation 
mode. The case study indicates the oscillation frequency is 
insensitive to the damping of the system, and the formulation 
may also be applied to an inter-area model of a multi-machine 
power system. 
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